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Abstract. We investigate the way how the total mass sum of neutrinos can be constrained from 
the neutrinoless double beta decay and cosmological probes with cosmic microwave background 
(WMAP 3-year results), large scale structures including 2dFGRS and SDSS data sets. First we 
discuss, in brief, on the current status of neutrino mass bounds from neutrino beta decays and cosmic 
constrain within the flat ACMD model. In addition, we explore the interacting neutrino dark-energy 
model, where the evolution of neutrino masses is determined by quintessence scalar filed, which is 
responsable for cosmic acceleration today. Assuming the flatness of the universe, the constraint we 
can derive from the current observation is £ m v < 0.87eV at the 95 % confidence level, which is 
consistent with £ m v < 0.68eV in the flat ACDM model. Finally we discuss the future prospect of 
the neutrino mass bound with weak-lensing effects. 
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NEUTRINO MASS BOUNDS FROM THE Ov/3/3 DECAYS AND 
COSMOLOGICAL PROBES WITH ACMD MODEL 

The existence of the tiny neutrino masses qualifies as the first evidence of new physics 
beyond the Standard Model. The answers to the hot questions on (1) whether neutrinos 
are Dirac or Majorana fermions ?, (2) the mass hierarchy pattern (normal or inverted 
hierarchy type ?), (3) the absolute value of the neutrino mass, will provide us the 
additional knowledge about the precise nature of this new physics, have the potential 
to unravel some of the deepest and most long-standing mysteries of cosmology and 
astrophysics, such as the origin of matter, the origin of heavy elements, and even the 
nature of dark-energy. 

The golbal analysis of the solar and KamLAND data [1] and super- Kamiokande 
atmospheric data [2] provide the two independent neutrino mass-squared differences: 
Am 2 sol = Am 2 = (7.9j£|) x l(T 5 eV 2 , \Am 2 atm \ = \Am 2 a \ = (2.6 ±0.2) x l(T 3 eV 2 , and 
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mixing angles: 9 l2 = 3 = (33.7 ± 1.3)°, 2 3 = 0i = (43.3+^)° and 13 = 2 < 5.2°. 
Above results tell us a substential eviden ce that t he three known neutrinos have a com- 
bined mass (E = L^i m v,i) of at least ^\Am a \ 2 ~ 0.05 eV. Since neutrino oscillations 
are only sensitive to mass-squared differences, three possible neutrino mass spectrum 
are allowed as: mi C m2 < /«3 (normal hierarchy), orm3<miCm2 (inverted hierar- 
chy), or mi ~ m 2 ~ m 3 (quasi-degenerate), depending on whether the light eigenmass is 
close to or 3> •y/jArafj", respectively. 

The nature of spectrum is important to neutrino mass model-building, the combination 
of neutrinos to dark-matter, and the viability of observing neutrinoless double beta-decay 
(Ov/3/3) if neutrinos are Majorana [3, 4, 5]. 

There are three well known ways to get the direct information on the absolute mass 
of neutrinos by using: Tritium /3 -decay experiment, neutrinoless double beta decay 
experiment, and astrophysical observations. 



(A) Neutrinoless Double Beta Decays: 

The standard method for the measurement of the absolute value of the neutrino mass 
is based on the detailed investigation of the high-energy part of the /3 -spectrum of the 
decay of tritium: 

3 H — > 3 He + e~ + v e (1) 

This decay has a small energy release (Eq ~ !S.6keV) and a convenient life time {T X j 2 = 
12.3 years). Since the flavour eigenstates are different from mass eigenstates in neutrino 
sector, in general, electron neutrino can be expressed as 

VeL = Y^UeiViL, (2) 

where v, is the field of neutrino with mass m,, and U is the unitary mixing matrix. 
Neglecting the recoil of the final nucleus, the spectrum of the electrons is given: 

and the resulting spectrum can be analyzed in term of a single mean-squared electron 
neutrino mass 

(m p ) 2 = Y,mj\U ej \ 2 = mj\U el \ 2 + m 2 2 \U e2 \ 2 + m 2 3 \U e3 \ 2 (4) 
j 

If the neutrino mass spectrum is practically degenerate: mi ~ m 2 ~ m 3 , the neutrino 
mass can be measured in these experiments. Present-day tritium experiments Mainz[6] 
and Troitsk[7] gave the following results: 

m\ = (-1.2±2.2±2.1)eV 2 (Mainz), (5) 
= (-2.3±2.5±2.0)eV 2 (Troitsk). (6) 



This value corresponds to the upper bound 



mi < 2.2eV (95%C.L.) (7) 

Another useful method is by using the neutrinoless double beta decay. The search for 
neutrinoless double /3 -decay 

(A,Z)^(A,Z2) + e-+e- (8) 

for some even-even nuclei is the most sensitive and direct way of investigating the 
nature of neutrinos with definite masses. In this process, total lepton number is violated 
(AL = 2) and is allowed only if massive neutrinos are Majorana particles. The rate of 
Ov/3/3 is approximately 

i = Gov(£/3/3,Z) l M «v| 2 (m^) 2 , (9) 

A l/2 

where G 0v is the phase space factor for the emission of the two electrons, Mq v is nuclear 
matrix elements, and < mpp > is the effective majorana mass of the electron neutrino: 

(^) = \%Uem\ (10) 

We can write eq.(10), for normal and inverted hierarchy respectively, in terms of mixing 
angles and A 2 = m\-m\ = (7.9±||) • 1(T 5 eV 2 , A a = ±(m 2 - m 2 ) ~ ±(2.6 ± 0.2) • 

\0~ 3 eV 2 and CP phases as follows: 



(m ee ) = 
(m ee ) = 



c\c\m x + c 2 2S 2 3 e^^A 2 + m 2 + S 2 e^^A 2 + m 2 



s\m x + c 2 S 2 e^^A 2 -A 2 + m 2 +c 2 S 2 e^ 



, ( normal hy erarchy ) ; 
(inverted hyerarchy ) .(11) 



From above relations, we can have the correlation plot between m/ ;g ^ and \tnpp\ with 
current observed data sets of mixing angles and A 2 a from neutrino oscillation experi- 
ments. However, Ov/3/3 decays have not yet been seen experimentally. 

The most stringent lower bounds for the time of life of Ov/3/3 -decay were obtained in 
the Heidelberg-Moscow[8] and IGEX[9] 16 Ge experiments: 



T^ 2 > 1.9 • 10 25 years (90%C.L.) Heidelberg - Moscow, (12) 

r 0v 

1/2 



T,°r 2 > 1.57 -\0 25 years (90%C.L.) IGEX. (13) 



Taking into account different calculation of the nuclear matrix elements, from these 
results the following upper bounds were obtained for the effective Majorana mass: 



\m 



00 1 < (0.35 - 1.24) eV (14) 



TABLE 1. The current upper limits on effective Majorana neutrino mass 
\mpp \ and the sensitivities of the future Ovj3/3-decay experiments. We used 
the matrix elements M 0v with reduced uncertainty [12]. T®J 2 denotes the 
current lower limit on the Ov/3/3 -decay half-life or the sensitivity of planned 
0vj3 j3 -decay experiments. 



Nucl. M 0v 


T®J 2 (years) Experiment \mpp (eV) 


16 Ge 2.40 


1.9 • 10 25 Hiedelberg-Moscow 0.55 
3 • 10 27 Majorana 0.044 
7 • 10 27 GEM 0.028 
1 • 10 28 GENIUS 0.023 


WQ Mo 1.16 


6.0 • 10 22 NEM03 7.8 
4 • 10 24 NEM03 0.92 
1 • 10 27 MOON 0.058 


m Te 1.50 


1.4 10 23 CUORE 3.9 
2 10 26 CUORE 0.10 


l36 Xe 0.98 


1.2 • 10 24 DAMA 2.3 

3 • 10 26 XMASS 0.10 
2 • 10 27 EXO(lt) 0.055 

4 • 10 28 EXO(lOt) 0.012 



Many new experiments (including CAMEO, CUORE, COBRA, EXO, GENIUS, MA- 
JORANA, MOON and XMASS experiments) on the search for the neutrinoless double 
/3 -decay are in preparation at present. In these experiments the sensitivities 

\m p p\ ~ (0.1-0.015)eV (15) 

are expected to be achieved. The detail upper limit of |m««| and the sensitivities of the 
future Ov/3/3 -decay experiments are summerized in table 1. It is very difficult to confirm 
the normal hierarchy pattern of neutrino mass when m\ < 1.7 • !0~ 3 eV, however for the 
inverted case, it can be detected if mj, < 8.9 • 10~ 3 eV and m ee > O.OlleV . 



(B) Cosmological Constrains within the Standard Cosmology: 

Within the standard cosmological model, the relic abundance of neutrinos at present 
epoch was come out straightforwardly from the fact that they follow the Fermi-Dirac 
distribution after freeze out, and their temperature is related to the CMB radiation 
temperature Tqmb today by T v = (4/1 \) 1 ^Tcmb with Tqmb = 2.726 K, providing 

„ - 6 ^ 3 ) T 3 nt* 

v — i i _2 1 CMBi V 1 ") 

where £(3) ~ 1.202, which gives n v ~ 1 12cm -3 for each family of neutrinos at present. 
By now the massive neutrinos become non-relativistic, and their contribution to the mass 



density (Q. v ) of the universe can be expressed as 



Q v h 2 = 



E 



93.14eV 



(17) 



where E stands for the sum of the neutrino masses. In this relation, we include the effect 
of three neutrino oscillation [10]. We should notice that when obtaining the limit of 
neutrino masses one usually assumes: 

• the standard spatially flat ACDM model with adiabatic primordial perturbations, 

• they have no non-standard interactions, 

• neutrinos decoupled from the thermal background at the temperatures of order 1 



These simple conditions can be modified from several effects: due to a sizable neutrino- 
antineutrino asymmetry, due to additional light scalar field coupled with neutrinos [13], 
and due to the light sterile neutrino [14]. However, analysis of WMAP and 2dFGRS data 
gave independent evidence for small lepton asymmetries [15, 16], and such a scenario 
with a light scalar field is strongly disfavored by the current CMB power spectrum data 
[17]. We will not therefore take into account such non-standard couplings of neutrinos 
in the following. In addition, current cosmological observations are sensitive to neutrino 
masses O.leV < E < 2.0eV. In this mass scale, the mass-square differences are small 
enough and all three active neutrinos are nearly degenerate in mass. Therefore we 
take the assumption of degenerate mass hierarchy. Even if we consider different mass 
hierarchy pattern, it will be very difficult to distinguish such hierarchy patterns from 
cosmological data alone [18]. 

After neutrinos decoupled from the thermal background, they stream freely and their 
density perturbations are damped on scale smaller than their free streaming scale. Conse- 
quently the perturbations of cold dark matter (CDM) and baryons grow more slowly be- 
cause of the missing gravitational contribution from neutrinos. The free streaming scale 
of relativistic neutrinos grows with the hubble horizon. When the neutrinos become non- 
relativistic, their freestreaming scale shrinks, and they fall back into the potential wells. 
The neutrino density perturbation with scales larger than the freestreaming scale resumes 
to trace those of the other species. Thus the free streaming effect suppresses the power 
spectrum on scales smaller than the horizon when the neutrinos become non-relativistic. 
The co-moving wavenumber corresponding to this scale is given by 



for degenerated neutrinos, with almost same mass m v . The growth of fourier modes with 
k > k nr will be suppressed because of neutrino free-streaming. The power spectrum of 
matter fluctuations can be written as 



MeV. 




(18) 



P m (k,z)=P*(k) T 2 (k,z), 



(19) 



where P*(k) is the primordial spectrum of matter fluctuations, to be a simple power law 
P*(k) =Ak n , where A is the amplitude and n is the spectral index. Here the transfer 



TABLE 2. Recent cosmological neutrino mass bounds (95% C.L.) 



Cosmological Data Set II E bound (2d) II References 



CMB (WMAP-3 year alone) 


< 2.0 eV 


Fukugita et al.[21] 


LSS[2dFGRS] 


< 1.8 eV 


Elgaroy et al. [22] 


CMB + LSS[2dFGRS] 


< 1.2 eV 


Sanchez et al.[23] 


tt 


< 1.0 eV 


Hannestad[24] 


CMB + LSS + SNla 


< 0.75 eV 


Barger et al.[25] 


tt 


< 0.68 eV 


Spergel et al.[26] 


CMB + LSS + SNla + BAO 


< 0.62 eV 


Goobaret al.[27] 


tt 


< 0.58 eV 




CMB + LSS + SNla + Ly-a 


< 0.21 eV 


Seljaket al.[28] 


CMB + LSS + SNla + BAO + Ly-a 


< 0.17 eV 


Seljaket al.[28] 



function T(k,z) represents the evolution of perturbation relative to the largest scale. If 
some fraction of the matter density (e.g., neutrinos or dark energy) is unable to cluster, 
the speed of growth of perturbation becomes slower. Because the contribution to the 
fraction of matter density from neutrinos is propotional to their masses (Eq. (17)), the 
larger mass leads to the smaller growth of perturbation. The suppression of the power 
spectrum on small scales is roughly proportional to f v [19]: 

(20) 

where f v = Cl v / Q.m is the fractional contribution of neutrinos to the total matter density. 
This result can be understood qualitatively from the fact that only a fraction (1 — f v ) of 
the matter can cluster when massive neutrinos are present [20]. Analyses of CMB data 
are not sensitive to neutrino masses if neutrinos behave as massless particles at the epoch 
of last scattering. According to the analytic consideration in [29], since the redshift when 
neutrino becomes non-relativistic is given by 1 +z nr = 6.24 • l0 4 Q v h 2 and z rec = 1088, 
neutrinos become non-relativistic before the last scattering when £l nu h 2 > 0.017 (i.e. 
E > \ .6eV). Therefore the dependence of the position of the first peak and the height of 
the first peak on Q. v h 2 has a turning point at Q. v h 2 ~ 0.017. This value also affects CMB 
anisotropy via the modification of the integrated Sachs-Wolfe effect due to the massive 
neutrinos. However an important role of CMB data is to constrain other parameters that 
are degenerate with E. Also, since there is a range of scales common to the CMB and 
LSS experiments, CMB data provides an important constraint on the bias parameters. 
We summarize some of the recent cosmological neutrino mass bounds within the flat- 
ACDM model in table 2. 



NEUTRINO MASS BOUNDS IN INTERACTING 
NEUTRINO-DARK ENERGY MODEL 

With our previous works [30, 31, 32], we investigate the cosmological implication of an 
idea of the dark-energy interacting with neutrinos [33, 34]. For simplicity, we consider 
the case that dark-energy and neutrinos are coupled such that the mass of the neutrinos 



is a function of the scalar field which drives the late time accelerated expansion of the 
universe. 

In our scenario, Equations for quintessence scalar field are given by 



$ + 2^+a 2 ^%M = 0, (21) 
VM) = (22) 

m) = a_4 /(03^ 2+a2m v^)^)' (23) 

mv (0) = m ; /4 , (24) 

where V (0) is the potential of quintessence scalar field, Vi((j>) is additional potential due 
to the coupling to neutrino particles [34, 35], and m v (<j>) is the mass of neutrino coupled 
to the scalar field, where we assume the exponential coupling with a coupling parameter 
p\ is |, where the dot represents the derivative with respect to the conformal time x. 

Energy densities of mass varying neutrino (MaVaNs) and quintessence scalar field are 
described as 

Pv = a- 4 J ■^J^\/q 2 + a 2 m 2 v f (q) , (25) 

/ (t% ,l \ 2 /o ^ ' (26) 

f + V(<i>), (27) 

f-V($). (28) 

From equations (25) and (26), the equation of motion for the background energy density 
of neutrinos is given by 

p v + 3Jf(p v +P V ) = ^^HPv ~ 3A0 • (29) 

Here we consider three different types of the quintessence potential V(0): (1) inverse 
power law potentials (Model I), (2) SUGRA type potential models (Model II), (3) 
exponential type potentials (Model III), which are given by, respectively: 

V{4>)=M*(^y ; M 4 (^y^ 2 / 2M P> ; mV"^. (30) 

The coupling between cosmological neutrinos and dark energy quintessence could mod- 
ify the CMB and matter power spectra significantly. It is therefore possible and also 
important to put constraints on coupling parameters from current observations. For this 
purpose, we use the WMAP3 [36, 37] and 2dFGRS [38] data sets. 
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FIGURE 1. (Left panel):Contours of constant relative probabilities in two dimensional parameter 
planes for inverse power law models. Lines correspond to 68% and 95.4% confidence limits; (Right 
panel):Same as Fig.l-a, but for exponential type models. 



The flux power spectrum of the Lyman-a forest can be used to measure the matter 
power spectrum at small scales around z < 3 [39, 40]. It has been shown, however, that 
the resultant constraint on neutrino mass can vary significantly from £m v < 0.2eV to 
0.4eV depending on the specific Lyman-a analysis used [41]. The complication arises 
because the result suffers from the systematic uncertainty regarding to the model for 
the intergalactic physical effects, i.e., damping wings, ionizing radiation fluctuations, 
galactic winds, and so on [42]. Therefore, we conservatively omit the Lyman-a forest 
data from our current analysis. 

Because there are many other cosmological parameters than the MaVaNu parameters, 
we follow the Markov Chain Monte Carlo(MCMC) global fit approach [43] to explore 
the likelihood space and marginalize over the nuisance parameters to obtain the con- 
straint on parameters we are interested in. Our parameter space consists of 

P={Qtf?,ClJ?-,H,T;,A s ,n s ,mi,a,P) , (31) 

where co^h 2 and Q. c h 2 are the baryon and CDM densities in units of critical density, H is 
the hubble parameter, T is the optical depth of Compton scattering to the last scattering 
surface, A s and n s are the amplitude and spectral index of primordial density fluctuations, 
and (m,, a, /3) are the parameters of MaVaNs. As an example, allowed parameter's space 
are shown in Figs.(l) for the model I and III. In these figures we do not observe the strong 
degeneracy between the introduced parameters. This is why one can put tight constraints 
on MaVaNs parameters from observations. For both models we consider, larger a leads 
larger w at present. Therefore large a is not allowed due to the same reason that larger w 
is not allowed from the current observations. We find no observational signature which 
favors the coupling between MaVaNs and quintessence scalar field, and obtain the upper 
limit on the coupling parameter as shown in table 3. 



P < 0.46,0.47,0.58 (la); [1.12, 1.36, 1.53 (2a)], 



(32) 



TABLE 3. Global analysis data within 2a deviation for different types of the quintessence 
potential. 



Quantites 


Model I 


Model II 


Model III 


WMAP-3 data (ACDM) 


a 


<4.38 


0.10- 11.82 


< 1.41 


— 


P 


< 1.12 


< 1.36 


< 1.53 


— 




2 09-2 36 


2.09-2.35 


2 08-2 34 


2 23 + 07 


0. CDM h 2 [\Q 2 ] 


9.87-12.30 


9.85-12.40 


9.84-12.33 


12.8±0.8 


Ho 


58.39-72.10 


58.55-71.70 


58.99-71.58 


72 ±8 


7 

t-'re 


6.13-14.94 


4.00-14.78 


6.64-14.78 




n s 


0.92-0.99 


0.92-0.98 


0.92-0.98 


0.958 ±0.016 


A s [10 10 ] 


18.25-23.41 


18.20-23.32 


18.33-23.27 




il e [10 2 ] 


57.43-75.60 


57.59-75.02 


58.45-75.05 


71.6±5.5 


Age 1 Gyrs 


13.59-14.40 


13.59-14.35 


13.61-14.36 


13.73 ±0.16 


rw/z 2 [io 2 ] 


<0.95 


<0.91 


<0.84 


< 1.97(95%CX.) 




0.031-0.143 


0.028-0.139 


0.032-0.140 


0.089 ±0.030 
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FIGURE 2. Examples of the total mass contributions in the matter power spectrum in Model I (Left 
panel) and Model III (Right panel). For both panels we plot the best fitting lines (green dashed), lines 
with larger neutrino masses M v = 0.3 eV (blue dotted) and M v = 1 .0 eV (cyan dot-dashed) with the other 
parameters fixed to the best fitting values. Note that while lines with M v — 0.3 eV can fit to the data well 
by arranging the other cosmological parameters, lines with M v — 1 .0 eV can not. 



and the present mass of neutrinos is also limited to 

^v/^oday < 0.0044, 0.0048, 0.0048 (la); [0.0095, 0.0090, 0.0084 (2 a)], (33) 

for models I, II and III, respectively. When we apply the relation between the total sum 
of the neutrino masses M v and their contributions to the energy density of the universe: 
Q v /z 2 = M v /(93.14eV), we obtain the constraint on the total neutrino mass: M v < 
0.45 eV (68.5%C.L.) [0.87 eV (95%C.L.)] in the neutrino probe dark-energy model. The 
total neutrino mass contributions in the power spectrum is shown in Fig 2, where we can 
see the significant deviation from observation data in the case of large neutrino masses. 
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FIGURE 3. Forecast 2a sensitivities to the total neutrino mass from future cosmological experiments, 
compared to the values in agreement with present neutrino oscillation data. Left panel: sensitivities ex- 
pected for future CMB experiments (without lensing extraction), alone and combined with the completed 
SDSS galaxy redshift survey. Right panel: sensitivities expected for future CMB experiments including 
lensing information, alone and combined with future cosmic shear surveys. Here CMBpol refers to a 
hypothetical CMB experiment roughly corresponding to the Inflation Probe mission. (Fig. from Ref.[46]). 



Beyond the scope of our current analysis, there are other possibilities in cosmological 
probes of neutrino masses: 

• the evolution of cluster abundance with redshit may provide further constraints on 
neutrino masses, 

• the Lyman- a forest provides constaints on the matter power spectrum on scale of 
k ~ 1 hMpc~ l , where the effect of massive neutrinos is most viable, 

• Deep and wide weak lensing survey will make it possible, in the future, to perform 
weak lensing tomography of the matter density field. 

As shown in Fig. 3, the combination of weak lensing tomography and high-precision 
CMB -polarization experiments may reach sensitivities down to the lower bound of 0.06 
eV on the sum of the neutrino masses [44, 45, 46]. In this case, normal hierarchy pattern 
will be detectable. 
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